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CORRECTIONS OF 0{a 2 s ) TO THE DECAY OF AN INTERMEDIATE-MASS 

HIGGS BOSON INTO TWO PHOTONS* 

MATTHIAS STEINHAUSER 



The QCD correction of 0(a 2 s ) to the decay of the Standard Model Higgs boson 
into two photons is presented. We consider the contribution coming from dia- 
grams with a heavy top quark as virtual particle. The first three terms of the 
expansion in the inverse top mass is calculated. Expressing the result through 
the on-shell top mass M t , we find large coefficients for the power-suppressed 
terms whereas in the MS scheme the coefficients are tiny. 



The Higgs boson is the only undetected particle in the Standard Model (SM). Now 
that top quark is well established the discovery of the Higgs boson is one of the most 
important aims in high-energy physics. Current experiments at LEP ruled out a Higgs 
boson with mass Mh < 65.6 GeV □ via Bjorken's process e + e~ — ► Z — > ffH. The next 
generation of colliders will improve this limit. Despite the fact that for an intermediate- 
mass Higgs boson with M\y < Mh < 2My/ the decay into a b quark pair will be the 
dominant process, a detection in this channel is unlikely with the LHC because of the large 
background. A very promising channel for the discovery especially for Mh < 130 GeV is 
the decay into two photons, though the branching ratio is of (9(10 -3 ). T(H — > 77) also 
enters the cross section for Higgs-boson production through photon-photon fusion, which 
is the dominant process in future high energy e + e~ colliders. Therefore it is important 
to calculate higher order corrections to this process. In this note we consider the 0(a 2 s ) 
corrections to T(H — > 77) coming from a top quark which is heavy compared to the 
Higgs boson. Besides the leading term which is independent of the top-quark mass, we 
also calculate the first two power-suppressed corrections. Throughout this paper we work 
in the framework of dimensional regularisation with D = 4 — 2e. Large intermediate terms 
are treated with the help of FORM i. 

The coupling ^77 does not exist at Born level, so the decay of the Higgs boson into 
photons is a loop-induced process. At the one-loop level either fermions or W bosons may 
be present in the loop, and the decay rate reads 
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with Tf = M H /4m 2 and Tw = M\jAMy^. In the one-loop order Aw is dominating the 
correction. The contribution from light fermions is suppressed by their mass. We will 
focus our attention to A t {r t ) and consider QCD corrections of 0(a 2 ) to this amplitude. 
It is convenient to write 



A t (T t ) = A^ + ^ + m 2 AV + ..., (2) 



7T \ 7T 



where a s = a^(fx) throughout this paper. 

There is no decoupling in the limit M t — > 00 because the Hit coupling is proportional 
to the top quark mass. The exact one-loop result is well-known and reads for r t < 1 ii 



4 0) = A 



k V 1 + V 1 - 7) arcsin2 ^) 



j / 7 2 2 26 o 512 4 1216 «, /nl fi A 

= if 1 + — r 4 + — r 2 + r, 3 + r 4 + r, 5 + O (rf) (3) 

V 30 21 * 525 * 17325 * 63063 * Kt 7 V ; 

with A 4 = N c ^Ql and t> = 2~ 1//4 G F 1//2 . The second line of Eq. (|3|) contains the expansion 
for large M t . 

At the two-loop level also the QCD corrections are at hand. Whereas for arbitrary 
values of Mh and M t only a numerical result is available i in the limit of a heavy top 
quark also analytical results exist i. We repeated this calculation and even extended it 
to get the first six terms in the expansion for a heavy quark. The result is given by 
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In order to get an impression of the quality of the expansion for large M t we plot in 
Fig. [I] the one- and two-loop results for A t against r t . The dotted lines are the result for 
infinitely large top mass. In the other curves successively higher orders of r t are included. 
The two vertical lines indicate M H = M w and M H = 2M W , respectively. Also the exact 
results are shown (full line). At two loops the result in the on-shell scheme (dashed line) 
differs from the one in the MS scheme (dash-dotted line, fi 2 = rh 2 ; here the variable of 
the abscissa is f t (fi) = / '4rh 2 (fi)) . This numerical difference is compensated by an 
appropriate change of the one- loop result. It should be noted that in this context the 
transitition to the MS scheme is done by replacing the on-shell top mass M t through the 
running mass mt(/i). 

For Mh = V%M t ~ 250 GeV (i.e. Tt = 0.5) the curves including the r 4 5 term are 
practically indistinguishable from the exact result. But already the curves which only 
contain the r 2 term deliver a reasonable result and for Mh ~ 2M^ this correction is a 
very good approximation. For Mh ~ M\y the inclusion of the linear term only is enough 
to describe the exact result very well. We take these considerations as a motivation to 
calculate in the three-loop case the terms proportional to t 2 . 
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Figure 1: One- and two-loop corrections to the amplitude A t . Successively higher orders 
in r t are included. Dotted line: leading order result, dashed line: on-shell scheme; dash- 
dotted line: MS scheme, full line: exact result. The two vertical lines indicate Mh = My/ 
and Mh = 2Mw, respectively. 

There are different methods which allow one to calculate of T(H — ► 77) in the limit 
that a heavy virtual quark is present in the loop and the Higgs boson is in the intermediate 
mass range i. 

A very elegant method to get the leading order (i.e. the limit M t — > 00) result is 
provided by the use of a low-energy theorem (LET) for the Higgs boson SB!. The basic 
idea is that the coupling of the Higgs boson to other particles is proportional to the mass 
of these particles. Both for fermions and gauge bosons the substitution m ; — > m^l+H/v) 
generates the coupling of the particle with mass rrii to the Higgs boson. For a Higgs boson 
carrying zero four-momentum this is obtained by differentiating the respective two-point 
function with respect to the virtual-particle masses. To be more precise, in our case we 
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can write down the effective Lagrangian 



C 
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i + n^(o) 



(5) 



where F° 



is the field strength tensor and n°^(0) = Ify*(q 2 )/q 2 \ q=0 . il°^(g 2 ) is (the 
transversal part of) the photon two-point function containing the top quark as virtual 
particle. The superscript zero indicates that we are still dealing with bare quantities. 
Applying the LET we get an effective Lagrangian for the #77 vertex 
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(6) 



It is clear that after the differentiation the renormalization of the coupling constant a° s 
and of has to be performed. 

The method of the LET has been used very extensively in the past. Besides effective 
if 77 and H gg couplings also the Hbb, HWW and HZZ couplings where considered at 
the two- 0'i and even three-loop levels 00. In a recent work the theoretical background 
of the LET was considered 0. The big advantage of the LET is that much less diagrams 
have to be taken into account and the three-point amplitudes are simply generated by 
differentiating w.r.t the top-quark mass. 

The result for IT?^(0) is given by0 
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where N, 



c = 3, C F = 4/3, C A = 3, L = ln(/i 2 /(M°) 2 ), T = 1/2 and n f = 6 is the 
total number of quarks. If this expression is inserted into Eq. (^j) and the top-mass and 
coupling-constant renormalization is performed we receive the leading-order term for the 
amplitude of the Higgs-boson decay into a pair of photons 
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with L = ln(/x 2 /M 2 ). In this paper the diagrams with the Higgs boson and the photons 
attached to different fermion lines are not considered. They will be treated elsewhere S 
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For this reason in order to compare the result from the LET with the explicit calculation, 
the contributions from the double-bubble diagram coming from the differentiation w.r.t. 
the inner top mass have to be subtracted. This is already done in (^). 

A second method to compute the leading-order result is provided by the use of the 
so-called Fock-Schwinger gauge for the external photons. In coordinate space this reads 
as x^A^(x) = 0. It can be shown @ that from this condition the expansion of A^(x) in 
terms of local operators has the form 

M*) = ^^(0) + -Lx^(D ff F w (0)) + ... • (9) 

In the present case we are interested in attaching two photons to the diagrams pictured 
in Fig. |2|. If only the leading order is considered it is enough to keep the first term in 
Eq. (|]). Then the prescription derived in is as follows: (i) Route the momentum 
associated with the photon through the diagram to the Higgs boson; (ii) differentiate 
the affected propagators w.r.t. this momentum; (in) set the momentum equal to zero 
and (iv) proceed in the same way with the second photon. At the three- loop level there 
are eight tadpole diagrams which have to be considered. A small computer program 
written in FORM I applies the above algorithm and attaches two photons. Afterwards 
the integration is performed. With this method the same result as in Eq. (§) is obtained. 
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Figure 2: Tadpole diagrams relevant for the method based on the Fock-Schwinger gauge. 
The diagram with a second light quark loop and the one containing the ghost particle are 
not displayed. 

There is also a method which leads directly to the decay rate. Thereby one has to 
consider U H (q 2 ), the two point function for the Higgs boson. The imaginary part of 
U H (q 2 ) coming from a two photon cut has to be extracted. In order to get the width the 



5 



relation 

T(H - 77) = -^lmU H (M 2 H ) (10) 

can be used. In our case, it is convenient to consider double-triangle diagrams where 
in both triangles top quarks (and possibly gluons) are present. In the limit that the 
external momentum q 2 = Mfj is less than 4M 2 , the only imaginary part comes from the 
two-photon cut. With this procedure we where able to reproduce the results of Ref. i. 

The most obvious method is the expansion of the triangle diagrams in their exter- 
nal momenta q± and q 2 . One thus successively obtains higher-order terms in Mfj/M 2 
depending on how far the expansion is performed. The amplitude for the decay of the 
Higgs boson into two photons with polarization vectors e M (gi) and e u (q 2 ) has the following 
Lorentz structure: 

i i 

where has no contribution for on-shell photons. In Eq. (|TT|) the sum runs over all 
diagrams relevant for the decay H — > 77. Due to gauge invariance, we have J2i a t,i = 
— J2i It is easy to find projectors for a t ^ and b t} %. 

a t,i = 77^ — To (<?i<?2 9nv - giz^ - qx^v) , (12) 



>t.i 



p- 2) 

(D-2)( qi q 2 ) 



'-qiQ2 g^u + qivqip, + {D- l)g^g 2 v) • (13) 



In fact, we calculated both, a ti i and b ti , in order to have an additional check for the 
correctness of our result. A (2D dimensional) Taylor expansion of the triangle diagrams 
in the two external momenta q\ and q 2 leads to the following structure for a t7 i(qi, q 2 )'- 

00 

Mgi,g 2 )= £ J lm M)\qlT{M2) n . (14) 

i,m,n=0 

After projecting out the Cg 0n , we end up with massive three-loop diagrams with zero 
external momentum. For these integrals the technique is meanwhile well established Ellll. 
After the decomposition of the numerator in terms of the denominator, the scalar integrals 
are reduced via recurrence relations i to trivial ones or so-called master integrals, for 
which a hard three-loop calculation is necessary. 

At 0(a 2 ), altogether 2 x 79 diagrams (the "2" comes from the two possible orientations 
of the fermion line) have to be taken into account. In Fig. |3], three exemplary graphs are 
pictured. Setting Ca = 3, Cp = 4/3 and rif = 6, the result reads (A t = J2i a t,i) 



A, = A, 
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where the dots represent higher orders in t% and a s . The leading-order result of this 
method coincides with Eq. (|8]). To the knowledge of the author, this is the first time that 
the LET was tested at three-loop level. In Eq. (|15f) , the on-shell definition of the pole 
mass was used. Expressing the result in terms of the MS mass 0, fh t (fi), the result looks 
like: 
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with / = In fi 2 /mf(fi). The leading order terms coincide in both schemes. This is the 
consequence of the fact that in the limit of infinitely heavy top-quark mass A t does not 
depend on M t at all. In Fig. |, the results of 0(a 2 ) are plotted. The curves show a similar 
behaviour as in the one- and two-loop cases. 
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Figure 4: 0(a 2 ) corrections to H — > 77. The same notation is adopted as in Fig. [I]. 



From Eqs. (El) and (|T6|) the same observation can be made as for earlier calculations 
up to 0(a 2 ) 000, namely the coefficients of a s /n are much smaller in the MS scheme 



than in the on-shell scheme. Consider e.g., /x 



M 2 



and /i 



m 



2 ■ 



in Eqs. flU) and (0). 



Then, the ratio of the coefficients in front of (a s /7c) 2 T t and (a s /Tr) 2 T 2 is about 50 and 20, 
respectively. Another feature concerning the MS scheme is already visible in the two-loop 
result of Fig. [I]. There the expansion in 1/M t seems to converge more rapidly in the MS 
than in the on-shell scheme. Actually, the curves including corrections of order r t 3 , or 
r t 5 are practically the same as the one where only the 0(r 2 ) corrections are present. The 
dash-dotted curves in Fig. [| indicate a similar behaviour. 

To conclude, the 0(a 2 ) QCD corrections to the decay process H — > 77 for an in- 
termediate Higgs boson were calculated. We considered the contribution from a heavy 
top quark and evaluated the first three terms in the expansion in the inverse top mass. 
Inspired from the one- and two-loop case this should be at least up to M H rs 2M w an 
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excellent approximation to the exact result. The leading term was confirmed using a low 
energy theorem and a method based on the Fock-Schwinger gauge. 
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